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Abstract. Multi-TeV e+e" collisions provide with a large enough sample of if" 
bosons to enable measurements of its suppressed decays. Results of a detailed study of 
the determination of the muon Yukawa coupling at ^/s — 3 TeV, based on full detector 
simulation and event reconstruction, are presented. The muon Yukawa coupling can 
be determined with a relative accuracy of 0.04 to 0.08 for Higgs bosons masses from 
120 GeV to 150 GcV, with an integrated luminosity of 5 ab^^. The result is not 
affected by overlapping 77 — > hadrons background. 

PACS numbers: 13.66.Fg, 14.80.Bn 

Submitted to: J. Phys. G: Nucl. Phys. 
1. Introduction 

The detailed investigation of the Higgs sector is anticipated as one of the central themes 
of the accelerator particle physics program after the Higgs boson will have been observed, 
most likely at the LHC The study of the Higgs profile will have to determine whether 
the observed particle is indeed responsible for generating the mass of gauge bosons, 
quarks and charged leptons. A significant corpus of studies, outlining strategies to 
test the Higgs mechanism for gauge bosons and quarks with accuracies down to a few 
percent, has already been assembled [U [2] These rely on the combination of data from 
the LHC and an e~^e~ linear collider operating at centre-of-mass energies ^/s = 0.25 TeV 
-1.0 TeV. In this paper we discuss the feasibility of a test of the Higgs mechanism in 
the lepton sector, by verifying the scaling through a precise determination of the 
Qh^i^ coupling in the if" — > decay, in 3 TeV e^e~ collisions. At the LHC, the 

significance of an observation of the —>■ fi~^fi~ decay is estimated to be only ~ 3 a 
for 115 GeV < Mh < 130 GeV, even when combining the CMS and ATLAS data on 
both the gluon and weak boson fusion production channels for 300 fb^^ of integrated 
luminosity [3]. It further decreases for larger values of Mfj- e^e~ collisions below 
1.0 TeV should provide with determination of qhtt to a relative statistical accuracy 
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of ~ 0.035 for a 120 GeV Higgs boson but, at best, with only a signal and not with 
a precision measurement in the gH^J.^J. Yukawa coupling [5]. The WW fusion process, 
e~^e~ H^UgUe, offers a large Higgs cross section, provided high enough ^/s energies 
can be attained, to make this measurement possible. In this paper the detection of 
the — > fi'^fJ'~ decay in high energy e~^e~ collisions is studied using full Geant-4- 
based simulation and reconstruction. The paper is organised as follows. In section 2 we 
discuss the simulation of the e+e^ H'^Ue^e, /^^A*^ process and of the relevant 

physics and machine-induced background. Section 3 describes the detector simulation 
and section 4 has the results. 

2. e+e" ^ i/VeZ/e, /i+^u" at = 3 TeV 

Higgs boson production in e^e~ collisions proceeds through two main processes: the 
so-called Higgstrahlung, e^e~ H^Z'^, and the WW fusion, e~^e~ H^Vei'e reactions. 
While the associated production of the Higgs boson with a Z° has a number of 
advantages, including the possibility to tag the presence of the Higgs boson independent 
on its decay modes, WW fusion reaches the largest cross sections, due to its oc log 
rise, if large enough centre-of-mass energies can be obtained. This makes the WW 
fusion process of special interest for the study of rare Higgs decays, such as bb 
at large Higgs boson masses, and — >• /i"*"/!". But as the energy increases, the 
production cross section is more and more peaked along the beam axis. Above 
~ 3 TeV, the further increase in production cross section is effectively lost, due to the 
limited detector acceptance in the very forward region. The production cross section at 
3 TeV and 5 TeV is respectively 2.3 and 3.0 times larger compared to that at 1 TeV. 
However, restricting the acceptance to Higgs bosons produced at | cos6'| < 0.99 (0.92), 
the gains in cross section become 2.1 (1.7) and 2.4 (1.9), respectively. The possibility 
of achieving e~^e~ collisions at multi-TeV energies, introduced by the CLIC concept |3], 
offers a real opportunity for observing this decay and measuring gHfifj. with enough 
accuracy to perform a significant test of the scaling of the Yukawa couplings in the 
lepton sector over a broad range of Higgs boson masses [5] . Multi-TeV e~^e~ collisions 
become important not only for exploring the highest energy frontier but also for precision 
studies of suppressed Higgs decays. 

3. Data Analysis 

Signal events e+e^ H^Uei^e, —>■ ^'^ ^~ are generated with Pythia 6.205 [6], 
including beamsstrahlung effects, at several Higgs boson mass values from 120 GeV up 
to 160 GeV. Cross sections are computed using CompHep 4.4.0 |7] and the Higgs decay 
branching fractions with HDECAY At -/s = 3 TeV, the effective e~^e~ — > H'^Ue^e 
production cross section, accounting for initial state radiation, is 0.48 pb to 0.45 pb and 
the decay branching fraction, BR(if° — » fi^jj,^) = 2.56x10""^ to 6.47x10"^, for Mh in 
the range from 120 GeV to 150 GeV. The irreducible ^'v^vp background is generated 
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with CompHep 4.4.0 interfaced with Pythia 6.205. The total cross section of events 
with an invariant mass of the system, M^^, in the range 100 GeV < M^^ < 

400 GeV is 5.31 fb. We assume to operate CLIC at centre of mass energy of 3 TeV for a 
total integrated luminosity of 5 ab~^, which corresponds to ~ eight years (1 yr = 10^ s) 
of operation at its nominal luminosity of 6.5x10^^ cm~^ s~^. The CLIC anticipated 
parameters are given in [9], for this analysis we use the luminosity spectrum as in |10j . 
This provides with ~ 2.3x10^ e^e" —>■ H^Uei^e events, 620 to 150 of which are followed 
by the — > fi^fi~ decay, depending on the value of the Higgs bosons mass in the 
range from 120 GeV to 150 GeV. The machine induced background, which affects this 
analysis most, is due to parasitic 77 — > hadrons events occurring within the same time- 
stamp as a physics event. The anticipated average number of such events for the CLIC 
parameters is 3.2 BX~^, or 161 events within a 25 ns time-stamp. 77 collisions are 
simulated using the GUINEAPIG program pTl [T2] and the cross sections parametrised 
according to [13]. The hadronic final states are simulated using Pythia and overlayed 
to a signal event. The average visible energy deposited in the detector by 77 —>■ hadrons 
events is ~ 90 GeV/BX. 

Charged particles are reconstructed in a multi-layered main tracker, made of high 
resolution silicon strip detectors, and a precision vertex tracker in a 5 T solenoidal field. 
Discrete Si tracking, following the example of the CMS experiment at LHC, emerged as 
an appealing solution for CLIC, where local hit occupancy and two-track separation may 
limit the applicability of gaseous trackers p3l ITO] . A large magnetic field is beneficial 
for ensuring an excellent momentum resolution for high pt particle tracks, as well as for 
confining low momentum charged particles from machine-induced backgrounds to small 
radii. This analysis is based on an implementation of the SiD detector [15], developed 
as a concept for the International Linear Collider (ILC) project, for the reconstruction 
of 3 TeV e~^e~ events. The SiD concept has a five-layered Si main tracker which ensures 
a momentum resolution Spt/pf < 5xl0~^ and polar angle coverage down to 8.3 o. The 
following assumptions on the detector response are made: the detection efficiency is 
0.98 per layer, the point resolution 7 /im, which includes alignment effects, and the 
detector is capable of time-stamping hits with 25 ns resolution, corresponding to 50 
bunch crossings (BX). 

While the SiD is designed for the reconstruction of events at centre-of-mass energies 
below 1 TeV, its geometry is well suited for the study of simple events, as those of 
our signal, up to multi-TeV energies. Full event simulation is performed using the 
SiD model implemented in the Mokka 06-03 program p77| based on Geant-4 jl8j . 
Simulation results are saved in the Icio format [19j and used as input to the event 
reconstruction. This is performed using a set of dedicated processors implemented in 
the Marlin reconstruction and analysis framework [20j. The reconstruction starts from 
hits generated in the sensitive layers along the track. The track pattern recognition 
is performed from the outermost main tracker hits and extends inwards to the vertex 
tracker. After the first stage of pattern recognition hit bundles are fitted using a simple 
helix. A second stage pattern recognition is performed using hits left unassociated after 
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the first stage. Reconstructed particle tracks are requested to have more than 4 hits and 
associated hits on at least 70 % of the number of sensitive surfaces traversed, depending 
on their polar angle, 6, fit probability in excess of 0.01, | cos^| < 0.98, and impact 
parameter significance below 3.5 times the extrapolation resolution, in both coordinates. 
These cuts remove poorly measured particle tracks. The momentum resolution of the 
reconstructed muon tracks is Spt/p^ = (3.46 ± 0.10) xlO~^, in agreement with the 
SiD specifications. Events with two, oppositely charged muon candidate tracks with 
momentum exceeding 10 GeV and below 750 GeV and event energy in charged particles 
100 GeV < -Echargod < 1-5 TcV are selected. 

4. Results 

The invariant mass distribution of the selected di-muon pairs is shown in Figure [H The 
distribution is fitted with a Gaussian curve to describe the signal and an exponential 
background. The mass resolution aMn/^H for signal events is 0.0034, which corresponds 




Figure 1. Reconstructed invariant mass distribution for signal and background 
+ -^'missing cvents with Mii =130 GeV and 5 ab^^ of integrated luminosity. 

to 0.40 GeV for Mj:^=120 GeV. The signal region is defined by a ± 2 ctmh interval 
around the fitted mass peak position and a binned fit is performed to extract the 
number of signal and background events. Results are summarised in Table [1] We obtain 
'fff-T^T)^ = 0-086 and 0.160 for Mh = 120 GeV and 150 GeV, respectively. These 
results correspond to a determination of gHfifi with a relative accuracy of 0.04 to 0.08. A 
significant signal for the decay can be obtained up to Mh = 155 GeV, when its branching 
fraction is only 4x 10~^, due to the rapid turn on of the WW decay channel of the Higgs 
boson. The analysis is repeated by overlaying the 77 — > hadrons background on signal 
events, for Mh = 120 GeV/c^. No degradation of the reconstruction efficiency and 
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resolution is observed for the amount of hadronic background corresponding to 50 BX 
(see Figure [2]). 




Figure 2. Reconstructed invariant mass distribution for signal events with Mh 
=120 GeV without (histogram) and with (points with error bars) 77 — > hadrons 
background overlayed. 

5. Conclusions 

The determination of the muon Yukawa coupling gH/ifi in e~^e~ collisions at centre- 
of-mass energies of 3 TeV at CLIC has been studied through the e~^e~ — > H^Ue^^e, 
— > 11^ ii~ process using full simulation and event reconstruction. The process is 
observable for Higgs mass values up to 155 GeV and the gH^^^l Yukawa coupling can be 
determined with a relative statistical accuracy of 0.04 to 0.08 for masses from 120 GeV 
to 150 GeV with an integrated luminosity of 5 ab^^. The superposition of machine- 
induced 77 hadrons background corresponding to 50 bunch crossings overlayed on a 
single iX^n'w event does not affect these results. 
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